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fT^ ' Abstract. 

We start with a Riemann-Hilbert Problems (RHP) with canonical normalization whose sewing functions depends on two or 
. more additional variables. Using Zakharov-Shabat theorem we are able to construct a family of ordinary differential operators 

^SJ ' for which the solution of the RHP is a common fundamental analytic solution. This family of operators obviously commute 

provided their coefficients satisfy certain nonlinear evolution equations. Thus we are able to construct new classes of integrable 
nonlinear evolution equations. We illustrate the method with an example of a new type 4-wave interactions. Its Lax pair 
consists of operators which are both quadratic in the spectral parameter A and take values in the so{5) algebra. 
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^ 1. INTRODUCTION 

d , 

The N-wave interactions describe special types of wave-wave interactions ll32i [T], [TtI S 113. which play 
important role in various fields in physics. 

The topic quickly attracted mathematicians from spectral theory, dynamical systems. Lie algebras, Hamiltonian 
dynamics, differential geometry, see ll29l l26l l33l l34l [sl [TtI 01 and the numerous references therein. It attracted also 
^ , a number of physicists because they found important applications of these nonlinear evolution equations (NLEE) in 

' fluid mechanics, nonlinear optics, superconductivity, plasma physics etc. As a result many different approaches for 

investigating the soliton equations and constructing their Lax representations, soliton solutions, integrals of motion, 
Hamiltonian hierarchies etc. were developed, see psl 126111 111 [lllMIM III- 

In the present paper, which is a natural extension of [8], we propose an alternative approach to the same class of 
. equations using as a starting point the Riemann-Hilbert problem (RHP) I37l [3^[33l[34ll29l[36ll23ll : the importance of 

the canonical normalization of RHP was noticed in 1 11, 6]. Our aim is to show that this allows one to construct lings 
of commuting operators and in addition gives a tool to study their spectral propeities. 

In Section 2 below we start with some preliminaiies concerning the RHP and their reductions. In Section 3 we use 
jets of order 1 to reproduce well known results about the 4-wave equations, see Chapter 3 of Ii29ll . In Section 4 we use 
jets of order 2 which allows us to construct new types of integrable 4-wave interactions whose interaction terms contain 
', quadratic and cubic nonUnearities, as well as jc-derivatives. The last Section contains discussion and conclusions. 

.9. 

2. PRELIMINARIES 
2.1. RHP with canonical normalization 

Let us formulate the RHP: 

^ + (x,f,A) = ^"(jc,f,A)G(x,f,A), A^eK, lim ^ + (x,f,A) = l, (1) 

where ^*(jc,f , A) take values in the simple Lie group © with Lie algebra g. ^ + (jic,f , A) (resp. ^^(x,f , A)) is an analytic 
functions of A for ImA*^ > (resp. for ImA*^ < 0). For simplicity we consider particular type of dependence of the 
sewing function G{x,t,X) on the auxiliary variables: 

/^-A*^[y,G(x,f,A)]=0, i^~X''[K,G{^,t,X)]=Q. (2) 
ax at 



where A; > 1 is a fixed integer and J^K are linearly independent elements of the Cartan subalgebra 7, A' e f) C g. 
The canonical normalization of the RHP means that we can introduce the asymptotic expansion 

^±(;i;,f,A) =expe(x,f,A), Q{x,t,X) ^Y,Qk{x,t)X-''. (3) 

k=\ 

Since ^^(x,f , A) are group elements then all Qk{x^t) G 0. However, 

g(x,f,A) =^±(;c,f,A)7|±(x,f,A), 3C(;c,f,A) ^<^±(;c,f,A)/:|±(jc,f,A), (4) 
belong to the algebra g for any J and K from g. Since K also belongs to the Cartan subalgebra t), then 

[a(x,f,A),3C(x,f,A)] =0. (5) 

An important tool in our considerations plays the well known Zakharov-Shabat theorem iIstI [ssll formulated below 

Theorem 1. Let (^*(x,f, A) be solutions to the RHP (|7]) whose sewing function depends on the auxiliary variables x 
and t via eq. (O. Then ^^{x,t,X) are fundamental solutions of the following set of differential operators: 



+f/(x,f,A)^±(x,f,A) - A*[y,^±(x,f,A)] = 0, 

ox 

=!-|- +y(jc,f ,A)^±(x,f, A) - A^'[7i:,^±(x,f ,A)] 0. 
at 



(6) 



Proof. The proof follows the lines of 0371 13811 . We introduce the functions: 

g±(x,f,A) = /^|±(x,f,A) +A'^<^±(x,f,A)7|±(x,f,A), 

(7) 

/t±(x,f,A) = i-^^^{x,t,X)+X^^^{x,t,X)K^^{x,t,X), 
at 

where = ((^*)^', and using (|2]l prove that 

^^(x,f ,A) = g^(x,f, A), A;^(x,f,A) = A;^(x,f ,A), (8) 
which means that these functions are analytic functions of A in the whole complex A-plane. Next we find that: 

lim ^+(x,f,A) = A*y, \\mk^{x,t,X)^X''K. (9) 

and make use of Liouville theorem to get 

k 

1 k-l 



g+{x,t,X)^g (x,f,A) = A'^7-£f/z(x,f)A* 



(10) 

k 

k+{x,t,X) = k-{x,t,X) = X^K -Y,Vi{x,t)X''-' . 

1=1 

We shall see below that the coefficients Ui{x,t) and Vi{x,t) can be expressed in terms of the asymptotic coefficients Q 
in eq. ©. □ 

Lemma 1. The operators L and M commute: 

i^-i^ + [U{x,t,X)-X''J,V{x,t,X)^X''K]^0. (11) 
at ox 

k k 

Uix,t,X)^Y.^,{x,t)X''-', V{x,t,X) = Y,Vi{x,t)X''-'. (12) 

1=1 1=1 

Proof. The operators L and M ^ have common fundamental analytic solutions (FAS), i.e. they must commute. The 
eqs. ( fTTT i are an immediate consequence of (|6|l. □ 



where 



2.2. Jets of order k 



In what follows we will consider the jets of order A: of A) and ^{xjX), see (|5]). We introduce them by: 

a(;c,f,A)= (A'=^±(x,f,A)y,|±(x,f,A)) =A'^7-C/(;c,f,A), 
X{x,t,X) = i^X''£,^{x,t,X)Kl,'^{x,t,X)j =X''K-V{x,t,X). 

The subscript + used above means that we insert the asymptotic expansions of ^ * and their inverse (|3]l and cut off the 
terms with negative powers of X. 

Obviously U (x) E q can be expressed in terms of Q{x). In doing this we take into account (|5]) and obtain |[l9ll 

3(^,f,A)=/+£^ad^7, X{x,t,X)^K+'^^eid''QK, (14) 

k=i k=i 

and therefore for Ui we get: 

1 , 

Ui{x,t) = -adgj7, U2{x,t) = -adg^y- -adg|7 
U^{x,t) = -adg,/- ^ (adg^adg, +adg,adQ,)y- ^adg^y 

(15) 

t/i:(x,f) = -adg.y-^ £ adgadg/-^ £ adgade^ade,./ 

s+p=k s+p+r=k 

and similar expressions for V; (x, f ) with 7 replaced by K. 

2.3. Reductions of polynomial bundles 

An important tool to construct new integrable NLEE is based on Mikhailov's group of reductions ll27ll . Below we 
will use mainly Z2 and Za? with N > 2 reduction groups. The basic Z2-examples are as follows: 

AQ'{x,t,eX*)A = -Q{x,t,X), 
BQ*{x,t,eX*)B^Q{x,t,X), (16) 
CQ\x,t,-X)C ^ -Q{x,t,X), 

where £^ = 1 and A, B and C are elements of the group 65 such that A^ = B^ = = t. As, for the ZAf-reductions we 
may have: 

DB,^{x,t,coX)b = E,^{x,t,X), DQ{x,t,(oX)D = Q{x,t,X), (17) 
where o)" = 1 and D'^ = 1 . 

These relations allow us to introduce algebraic relations between the matrix elements of Q{x,t,X) which will be 
automatically compatible with the NLEE. The classes of inequivalent reductions of the A^-wave equations related to 
the low-rank simple Lie algebras are given in IA .10..12[[l3l[l3lfl4lfl6ll . 

3. THE CLASSICAL 4- WAVE INTERACTIONS 

The 4-wave interactions were discovered by Zakharov (see Chapter 3 of ll29ll ) who used a 4 x 4 Lax pair with an 
additional reduction which effectively reduced it to the subalgebra sp{4). It is well known that sp{4) is isomorphic to 
so{5), so it is a matter of taste which one to choose. We remind that the root system A of so{5) has 8 roots (see U^): 



a) A^'^'^x,t,eX*)A = ^^{x,t,X), 

b) B^+'*{x,t,eX*)B^^-{x,t,X), 

c) C^+'^(x,f,-A)C = |-(x,f,A), 



A = A+ U A_, A± = {±ai , ±a2, ±(«i + 0:2), ±(0:1 + 2a2)}, 



(18) 



where a\—e\— and — are the two simple roots. 

With this notations the Lax pair proposed by Zakharov takes the usual form for the A^-wave equations 

Li// = + (f/i(jc,f) + A7)v/(x,f,A) = 0, 
ax 

My/ = + (yi(x,f) + A/i:)v/(x,f,A) = 0, 
at 

where the potentials Ui and Vi are of the form: 

Ui{x,t) = [J,Qi{x,t)], Vi{x,t) = [K,Qi{x,t)], 

2l(x,f)= {qa{x,t)Ea+ Pa{x,t)E-a) 

We impose the natural reduction (fTSK ) which provide pa = £q*a- Thus the 4-wave equations take the form 



J. 



dQi 
dt 



K 



dQi 
dx 



[[J,QiUK,Qi{x,t)]]=0, 



or in components: 



2/(fli 




-2/(^1 -/72)^r^+2K-V'2M4 = 
OX 






0, 




dt 


— 2/^72^7— +2K'(v4M3 + M4V1) 

dx 


= 0, 




du4 
2ia 1 ^r— 
dt 


-.■7 '^"4 / \ 

— 2ibi— |-2k'(v2M3 — U2U1) 

dx 


= 0, 


2i{a\ 


du3 


— 2/(^1 + /72)^r— ~ 2ku2U4 — 
dx 


0, 



where mi = q'aj, M2 = g'az' "4 = ?ai+a2 ™d "3 = ?ai+2a2- 

4. NEW TYPES OF 4-WAVE INTERACTIONS - AN EXAMPLE 

Here we will give an example of a new type of integrable 4-wave equations. 
The Lax pair for these new equations will be provided by: 

= + {U2{x,t)+XUi{x,t) - X^J)\ir{x,t,X) = 0, 
dx 



.dxj/ 

'~dt 



My/ = + {V2{x,t)+XVi{x,t) - l^K)yr{x,t,l) ^ 0, 



where Uj(x,t) and Vj{x,t) are fast decaying smooth functions taking values in the Lie algebra so (5) 

Ui{x,t) = [J,Qi{x,t}], U2{x,t) = [J,Q2{x,t)] - ^ad^/, 
Vi{x,t) = [K,Qiix,t)], V2ix,t) - [K,Q2{x,t)]-]-^dl^K. 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



Here adg^X = [2i(-x,f)i^] ™d J and K are two linearly independent constant elements of the Cartan subalgebra of 
so{5). 

If we assume Qi (x,f ) and Q2{x,t) to be generic elements of so (5) then they can be parametrized by: 

Qlix,t)= £ iql,Ea+plE^a)+r\H,,+r\H,,^, 

e2(x,f)= {qiEa + piE-a) + r^iHe,+rlHe„ (25) 

aGA+ 

J = aiHei +a2He2 =diag(fli,fl2,0,-a2,-fli), K = b\He^ +b2He2 = Amg{bi,b2,Q,-b2,-b\), 



Next we impose on Qi {x,t) and Q2{x,t) the natural reduction (fT6h ) 

BoU{x,t,eV)''B(^^ =U{x,t,X), Bo = diag(l, £,!,£,!), e^ = l. 

As a result: 



(26) 



(27) 

which provide p^, = e(^^)*, pj^ = e(^a)*- Then eq. (|32] | will be a (rather complicated) system of 8 NLEE for the 8 
independent matrix elements ql^ and q^. 

However we can impose additional Z2 reduction condition 

D^^{x,t,-X)D^^^{x,t,X), DQ{x,t,-X)D = Q{x,t,X), 

where D = diag (1,— 1,1,— 1,1). More precisely this means: 



(28) 



Ql{x,t) = UiEe^-e2 +U2Ee2 + U^Ee^+e2 +Vl£'-e|- 



/ Ml M3 \ 

VI M2 M3 

V2 M2 

V3 V2 Ml 

V V3 v'l / 



/wi M4 \ 
' W2 ' 



(29) 



1^4 M4 

-W2 
V -V4 -wi ) 

J = aiHej +a2He2 =diag(ai,fl2,0,-fl2,-fli), K = b\He^ +b2He2 = diag (/7i,/72,0, -^2, -^1), 
Combining both reductions for the matrix elements of Qj{x,t) we have: 

vi = EMj, v2 = eM|, V3 = eM3, V4 = M4, 
The commutativity condition for the Lax pair (l23T l 

i(^+X^)^i(^+X^)+lU2 + XU,-X'j,V2 + XV,-X'K]^0 



dx 



dx 



dt J 



(30) 



(31) 



must hold identically with respect to X. One can check, that with the choice (1251 ) for Ui and V,, 1,2 the terms 
proportional to A^, and A^ vanish identically. The term proportional to A and the A-independent term vanish 
provided Q, satisfy the NLEE: 



dVi dUi , 

- + [U2,Vi] + [UuVi] = 0, 

dV2 dU2 , 

In components the corresponding NLEE: 

—2i{ai — 02)-^ + 2i{b\ — b2)—^ + JceMKeM^Ms — mim2 — 2M4) = 0, 

— 2;fl2^7— + 2ib2^r— — K'(M2e(|M3p — |mi P) +2M3M4 + 2£mTm4) = 0, 

at ax 

—2i{a\ +02)^7— +2i(b\ +^2)^7^ + K'M2(eM2M3 — M1M2 + 2uA — 0, 
at ax 

— 2ifli^r-^ +2ibi^r-^ + i— (— (2fl2 — fli)MiM2 + (2fl2 + fli)eM2M3) + ((2^2 — bi) ^ ^ 
at ax at ox 

-i{2b2 + bi)e^^-j^ - K-(2eM4(|Mi|^ - |M3p) + eMiM2(|Mip + 3|m3|^) - m3M2(3|miP + |m3|^)) =0. 



(32) 



(33) 



Let us now introduce 



As a result we get: 



U4 = U4 — ((fli — 02)^1 M2 + (fll +fl2)£M3M2)- (34) 



-2/(fli — ai)^ \-2i(bi — b2)^ U2{2a\U4 + efl2M2"3 + (2^1 ^ a2)u\U2) — 0, 

— 2/02^7— +2/^2^7— M2((2fli +a2)|M3 P ~ a2|Mi ~2k{u^Ua + £u*,U4 + uluXuj) — 0, 

at ox a\ 

-2/(fli +02)^7— +2/(^71 + /72)^r- H M2(e(2fli +fl2)M2"3 ~ fl2MiM2 + 2aiL'4) — 0, (35) 

31/4. dU4 iKdu\U2 iK£ du^u^ 

-2iai ^— + 2ibi ^— H 5 ^ 

at ax ai ox fli ox 

(2ef/4(|Ml| — |m3| ) + (eMiM2 — M3M2)((2fll ^ a2)lMl I + (2fll +fl2)|M3| )) = 0, 



(36) 



5. DISCUSSION AND CONCLUSIONS 

On the example of RHP for functions taking values in the 5(9(5)-group we have outlined the construction of families 
of commuting operators. Applied to jets of order 1 with this method reproduces the well known results for the 4-wave 
equations. 

Using jets of order 2 gives us the simplest nontrivial examples for new types of integrable 4-wave equation whose 
interaction terms contain quadratic and cubic nonlinearities, as well as x-derivatives. These equations also allow 
integrable extensions to three-dimensional space-time, see |8]. 

It is not difficult to obtain other new integrable 4-wave equations. To this end we can use: i) different power k of the 
polynomials U{x,t,X) and V{x,t,X) ii) different types of grading; and iii) different reductions of U and V. 

For example another integrable 4-wave interactions model can be obtained using Lax pair of degree 4 with respect 
to A: 

Lxj/ = + (U4(x,t) + XU3(x,t) + X^U2(x,t) + X^Ui(x,t) - X'^J)\l/(x,t,X) = 0, 
ox 

Mxff = + (V4(x,t) + XV^ix^t) + X *2V2(x,t) + XM{x,t) - X'^K)\(fixJ,X) ^ 0, 
ot 

where Uj{x,t) and Vj(x,t) for j = 1,2 are again as in ( l24] i: for j = 3 and j = 4 we have 

1 1 , 

U3{x,t) ^ [J,Q3{x,t)] - - (adgjadg, +adg2adQ|)7+ -adg^y, 

1 1 , 

Vs (jc, f ) = 63 f )] - 2 (^'l Gi ad 82 + ad 82 ad e, ) + - ad 

U4{x,t)^[J,Q4{x,t)]-^ £ ade,adQ/-i ^ adg.adg^.adQ^y - ^^ad^/, 

^(+7=4 ^ i+j+k=4 

U4(x,t)^[K,Q4(x,t)]~^ Y adQ_adQ^K-^ ^ adaadg^adg,/: - :^ad^^7i:, 

^i+;=4 ^i+j+k=4 

To keep the number of the independent functions down to 4 we will use in addition to (fT6h). an additional Z4- 
reduction of the form (fTTt which leads to the following parametrization of Qj: 

Ql{x,t) = UiEe^-e2 +U2Ee2 +^'l£'-c,+e2 +^'2£'-(!2i Qlixj) = M3£'t.,_e2 +^'3£'-C2' 

(3o) 

Q3[X,t) == U4Eej+e2+V4E-ei-e2, Q4-(x,t) == W4He^ +Z4He2, 

It jets U{x,t,X) can be view as elements of special co-adjoint orbits of the relevant Kac-Moody algebra, generated 
by the chosen grading of f. This allows one to define a hierarchy of Poisson brackets, see [25, 30, 31.1, which along 
with the conservation laws will provide the hierarchy of Hamiltonian structures of these NLEE. 



(37) 



Finally we can also apply Zakharov-Shabat dressing method iIstI [ssi l28i I21I1 for constructing their explicit (A^- 
soliton) solutions. Instead of solving the inverse scattering problem for L we would rather deal with a Riemann-Hilbert 
problem with canonical normalization. For polynomials of order k the contour on which the RHP is defined consists 
of k straight lines 4 : arg X = ni/ k passing through the origin. Of course, it may necessary to use dressing factors with 
more specific A -dependence. 

This approach can be used also to analyze the NLEE derived by Gel'fand-Dickey approach fslflsl]. It would provide 
the possibility to systematically construct the spectral decompositions that linearize the relevant NLEE ||7l M7D. Still 
more challenging is to study the soliton interactions of the new types of 4- and A^-wave equations. 
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